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Description 

[METHOD AND APPARATUS FOR 
ALGORITHM FUSION OF 

HIGH-RESOLUTION 
ELECTROCARDIOG RAMS] 

Background of Invention 

[0001] This invention relates to the analysis of electrocardio- 
grams (ECGs). More specifically, the invention relates to 
methods and apparatus for analyzing ECGs. 

[0002] | n hospitals or other health-care settings, it is frequently 
necessary to observe critical physiological conditions of a 
patient, including cardiovascular conditions. Such cardio- 
vascular condition data includes ECGs acquired using 
electrodes applied to the patient. 

[0003] a typical ECG management system comprises a database 
of ECGs plus applications software. The ECG management 
system receives ECG data from a multiplicity of instru- 
ments via a plurality of networks, analyzes that ECG data 



using various programs such as a serial comparison pro- 
gram, generates reports containing the results of the 
analysis, and then routes those reports to various systems 
and peripheral devices. 
[0004] The accuracy of any ECG analysis expert software tool is 
directly dependent upon the quality of the signal it ac- 
quires. In 1979, Marquette Medical Systems introduced an 
electrocardiograph that simultaneously acquired all of the 
leads from the 12-lead electrocardiogram. Prior to this 
time, all commercially available electrocardiographs could 
only acquire 3 leads at a time. Simultaneous recording 
was adopted so that the computer could use all signals 
from all 12 leads to properly detect and classify each QRS 
complex. The program also applied digital filters which 
removed power line noise and baseline sway. Typically, 
ECG data is acquired from the 12 leads for a period of 10 
seconds. 

[0005] All ECG analysis computer programs are composed of two 
parts: one which measures the waveforms, the other 
which does the interpretation based on these measure- 
ments. The main task of the measurement section is to 
find the location of the major reference points (that is, the 
onsets and offsets of the P, QRS and T complexes). After 



the onsets and offsets of the P, QRS, and T complexes 
have been demarcated, the waves within each complex are 
measured according to published standards. These ampli- 
tudes and durations result in a measurement matrix con- 
taining more than 1000 values. This is then passed to the 
criteria portion of the ECG analysis program so that it can 
generate an interpretation, including diagnostic state- 
ments referenced via a statement library. 
[0006] Recent progress in noninvasive electrocardiology has in- 
troduced more interesting features for use in diagnosis, 
including T wave alternans (TWA), QT dispersion (QTD), 
QT dynamicity (QTDN), short-time heart rate variability 
(HRV), etc. 

[0007] j wave AlternansAlternans is a subtle beat-to-beat 

change in the repeating pattern of an ECG that can be in- 
dicative of electrical instability of the heart and increased 
susceptibility to sudden cardiac death. Alternans results in 
an ABABAB . . . pattern of variation of the waveform shape 
between successive beats in an ECG waveform, and the 
level of cardiac stability is taken as a characterization of 
an individual's cardiac electrical stability. In humans, it has 
been found that the T wave is the best interval of the ECG 
complex for detecting alternans. That is, a level of varia- 



tion in the T waves of alternating beats is the best indica- 
tor of a patient's level of alternans. 

[0008] qt Dispersion The calculation of QT interval dispersion 

from the standard 12-lead ECC is a promising noninvasive 
measurement of inhomogeneity in myocardial repolariza- 
tion. The QT dispersion measurement is useful for differ- 
ent clinical applications, namely for assessment of pa- 
tients at risk for lethal ventricular arrhythmias and for dif- 
ferentiating patients with acute myocardial infarction 
(AMI) from those with noncardiac chest pain. To compute 
QT dispersion, one needs to measure QRS complex onset 
and T wave offset to calculate the QT interval for each 
lead. A computerized automatic method for calculating QT 
dispersion was presented by Xue and Reddy in "Algo- 
rithms for Computerized QT Analysis", Journal of Electro- 
cardiology, Vol. 30 Supplement, pp. 198-203 (1998). 

[0009] qt Dynamicitylt is known that prolongation of the QT in- 
terval may be a marker for sudden death. QT dynamicity is 
used to examine the beat-to-beat change of the QT inter- 
val due to changes in heart rate or T morphology. 

[0010] Heart Rate VariabilityThe autonomic nervous system is the 
part of the vertebrate nervous system that regulates in- 
voluntary action, such as the heart. A significant relation- 



ship exists between the autonomic nervous system and 
cardiovascular mortality, including sudden cardiac death 
(SCD). Heart rate variability (HRV) is a recognized quanti- 
tative marker of autonomic activity. Many commercial de- 
vices provide automated measurement of HRV. The term 
"heart rate variability" has become the conventionally ac- 
cepted term to describe variations of both instantaneous 
heart rate and RR intervals (i.e., the interval between con- 
secutive heart beats). Standards for the measurement of 
HRV are set forth in "Heart Rate Variability: Standards of 
Measurement, Physiological Interpretation, and Clinical 
Use", by the Task Force of the European Society of Cardi- 
ology and the North American Society of Pacing and Elec- 
trophysiology, Circulation, Vol. 93, No. 5, March 1, 1996, 
pp. 1043-1065. 
1 ] Late PotentialsLate potential is the small amplitude signal 
detected at the end of the QRS complex, so-called "late" 
potential. Since late potential has very low amplitude 
(usually lower than 50 uV), it is detected by using high- 
resolution ECG through a signal averaging technique. In 
this case "high resolution" means higher amplifier gain 
(around 1 uV/LSB vs. 5 uV/LSB (least significant bit)) and 
higher sampling frequency (1000 samples per second or 



more vs. 125-500 samples per second for conventional 
ECG). It has been proposed to combine T-wave alternans 
and late potentials for identifying high-risk patients with 
impaired left ventricular function. See Kondo et al., "Clini- 
cal usefulness of the combination of T wave alternans and 
late potentials for identifying high-risk patients with 
moderately or severely impaired left ventricular func- 
tions", Circulation, Vol. 65, No. 7, pp. 649-553 (2001). 
The late potentials were determined on signal-averaged 
ECGs. 

[0012] |ntra-QRS Analysislntra-QRS analysis examines small 

notches inside the QRS complex based on signal-averaged 
high-resolution ECGs. It has been shown that the signifi- 
cant notches found inside the QRS complex are correlated 
with ischemic events before and after PTCA. See Endt et 
al., "Identification of post-myocardial infarction patients 
with ventricular tachycardia by time-domain intra-QRS 
analysis of signal-averaged electrocardiogram and mag- 
netocardiogram", Med. Biol. Eng. Comput., Vol. 38, No. 6, 
pp. 659-665 (2000)Many other expert software tools are 
being used in cardiology. For example, high-resolution 
ECGs, also referred to as signal-averaged ECGs (SAECGs), 
have been widely used in cardiology clinical practice. Late 



potentials detection of SAECGs has been proven useful for 
predicting lethal ventricular tachyarrhythmia and ventricu- 
lar fibrillation (VT/VF), and P-wave SAECGs have been 
used for predicting atrial fibrillation (AF). During the pro- 
cess of high-resolution detection, 200-400 heart beats 
signals are acquired with higher sampling frequency 
(1000 Hz vs. 500 Hz for standard 12-lead ECG), and 
higher resolution (about 1 uV least significant bit vs. 5 uV 
for standard 12-lead ECG). In current cardiographs, the 
end results of high-resolution detection are averaged 
ECGs from three orthogonal leads: X, Y, Z, and some mea- 
surements like P wave and QRS duration from the aver- 
aged leads. In practice, also 12- lead electrodes are put 
on a patient"s body surface. 
[0013] with the 200-400 heart beats (5 to 10 minutes) of high- 
resolution signals from X, Y, Z leads plus standard 
12-lead ECGs, a cardiologist has much more information 
than the data acquired by regular 10-sec 12-lead ECGs. 
There is a need for methods of extracting more useful 
features from this increased volume of acquired data for 
more accurate and more efficient diagnosis of abnormal 

activity of the heart. 
Summary of Invention 



[0014] The present invention is directed to methods and appara- 
tus for multiple algorithm fusion of high-resolution ECG 
data. The disclosed embodiments use signals collected 
during high-resolution processing to not only detect new 
features such as late potential, T wave alternans (TWA), 
QT dispersion (QTD), QT dynamicity (QTDN), short-time 
heart rate variability (HRV), and intra QRS, but also to im- 
prove the quality of those features. For example, the re- 
producibility of QT dispersion measurements can be im- 
proved with many more heart beats of data. More impor- 
tantly, the combined features (feature fusion) from this 
group can improve the accuracy of the diagnosis in clini- 
cally relevant groups. For example, results from late po- 
tentials, T wave alternans and QT dynamicity analyses can 
be fused to predict ventricular tachyarrhythmia after my- 
ocardial infarction and atrial fibrillation; results from QT 
dispersion and intra-QRS analyses can be fused for is- 
chemia detection and prediction; and results from HRV, 
QT dispersion and QT dynamicity can be fused for pre- 
dicting sudden cardiac death after myocardial infarction. 

[0015] one aspect of the invention is a method for noninvasive 
ECG detection and diagnosis, comprising the following 
steps: acquiring high-resolution ECG data from a patient; 



processing the acquired data in accordance with two or 
more different ECG analysis algorithms; and deriving a 
prediction score for a particular clinical end point as a 
function of the respective results of those ECG analysis al- 
gorithms. 

[0016] Another aspect of the invention is a system for noninva- 
sive ECG detection and diagnosis, comprising: a multiplic- 
ity of electrodes applied to a patient, a data acquisition 
system for acquiring high-resolution ECG data from the 
patient; and a processor programmed to process the ac- 
quired data in accordance with two or more different ECG 
analysis algorithms, and then derive a prediction score for 
a particular clinical end point as a function of the respec- 
tive results of those ECG analysis algorithms. 

[0017] a further aspect of the invention is a method for noninva- 
sive ECG detection and diagnosis, comprising the follow- 
ing steps: (a) acquiring high-resolution ECG data from a 
patient using the X, Y and Z leads of a Frank lead system; 
(b) dividing the high-resolution ECG data into even beat 
and odd beat groups; (c) averaging beats in said even and 
odd beat groups separately; (d) determining the variance 
of T wave morphology of all even beats for each of said X, 
Y, Z leads; (e) determining the variance of T wave mor- 



phology of all odd beats for each of said X, Y, Z leads; (f) 
determining the variance of T wave morphology between 
even and odd averaged beats for each of the X, Y, Z leads; 
and (g) determining T wave alternans as a function of the 
results of steps (d)-(f) by analysis of variance. 
[0018] yet another aspect of the invention is a system for nonin- 
vasive ECG detection and diagnosis, comprising: a multi- 
plicity of electrodes applied to a patient, a data acquisition 
system for acquiring high-resolution ECG data from the 
patient using the X, Y and Z leads of a Frank lead system; 
and a processor programmed to process the acquired data 
in accordance with a T wave alternans algorithm compris- 
ing the following steps: (a) dividing the high-resolution 
ECG data into even beat and odd beat groups; (b) averag- 
ing beats in said even and odd beat groups separately; (c) 
determining the variance of T wave morphology of all 
even beats for each of said X, Y, Z leads; (d) determining 
the variance of T wave morphology of all odd beats for 
each of said X, Y, Z leads; (e) determining the variance of 
T wave morphology between even and odd averaged beats 
for each of the X, Y, Z leads; and (f) determining T wave 
alternans as a function of the results of steps (c)-(e) by 
analysis of variance. 



[0019] other aspects of the invention are disclosed and claimed 
below. 

Brief Description of Drawings 



[0020] FIG. 1 is a block diagram showing a system in accordance 

with one embodiment of the invention. 
[0021] FIG. 2 is a graph showing conventional ECG nomenclature. 

[0022] FIG. 3 is a drawing showing the standard 12-lead elec- 
trode system. 

[0023] FIG. 4 is a drawing showing the Frank lead electrode sys- 
tem. 

[0024] FIG. 5 is a flowchart of algorithm fusion for high- 
resolution ECG analysis in accordance with one aspect of 
the invention. 

[0025] FIG. 6 is a flowchart showing three-stage processing per- 
formed during high-resolution ECG analysis in accordance 
with another aspect of the invention. 

[0026] FIG. 7 is a diagram showing the application of a T wave 

alternans algorithm to high-resolution ECG in accordance 
a further aspect of the invention. 

[0027] Reference will now be made to the drawings in which sim- 
ilar elements in different drawings bear the same refer- 
ence numerals. 



Detailed Description 

[0028] Referring to FIG. 1, an electrocardiogram (ECG) system in 
accordance with one embodiment of the present invention 
comprises a set of electrodes 2 that are attached to the 
patient"s body, a data acquisition system 4 that receives 
and stores ECG waveform data output from the various 
electrodes 2, and a data processor 6 for processing the 
acquired ECG waveform data in accordance with ECG anal- 
ysis and algorithm fusion applications software. The data 
acquisition system 4 acquires samples of signals pro- 
duced by both orthogonal X, Y, Z leads and the standard 
12-lead system. A system operator can select the operat- 
ing mode of the processor via inputs to an operator inter- 
face 10. The results of the ECG analyses and subsequent 
algorithm fusion are sent from the processor to one or 
more output devices 8, which may include a display moni- 
tor, a printer, a storage medium, etc. 

[0029] The electrodes 2 are attached to a patient's skin and posi- 
tioned to detect electrical signals produced by the pa- 
tient's heart. The electrodes include an electrically con- 
ductive gel that contacts the patient's skin and conducts 
to the electrode any electrical signals that are present at 
the skin. The patient's heart produces an electrical signal 



that is referred to as an ECG waveform. The shape of a 
typical ECG waveform is shown in FIG. 2. The standard 
nomenclature identifies a P wave, a QRS complex, a T 
wave, and a U wave. The interval from the onset of the 
QRS complex to the end of the T wave is referred to as the 
QT interval. 

[0030] An ECG waveform for a single beat is typically referred to 
as a PQRST complex. The P wave appears at initiation of 
the beat and corresponds to activity in the atria, while the 
QRST complex follows the P wave and corresponds to 
ventricular activity. The QRS component represents the 
electrical activation of the ventricles, while the T wave 
represents the electrical recovery thereof. The ST segment 
is a relatively quiescent period. In humans, it has been 
found that the T wave is the best interval of the ECG com- 
plex for detecting alternans. That is, a level of variations 
in the T waves of alternating beats is a good indicator of a 
patient's cardiac electrical stability. 

[0031] The ECG signal produced by the patient's heart decreases 
as a function of the distance from the heart at which it is 
measured. Accordingly, an ECG signal detected by an 
electrode 2 will vary from the actual ECG signal based on 
the placement of the electrode relative to the heart. An 



accurate approximation of the actual ECG signal may be 
generated by combining signals from multiple electrodes 
having known placement relative to the heart. 

[0032] The ECG signal measured at the body surface may be rep- 
resented by modeling the heart as a three-dimensional 
dipole source that produces an electrical signal that varies 
based on the distance from the heart in the X, Y and Z di- 
rections. Use of the dipole vector model of the heart has 
lead to the development of clinical systems that measure 
the X, Y and Z components of the dipole source through a 
linear combination of several electrodes. One known XYZ 
system is the Frank lead system. 

[0033] | n accordance with one aspect of the present invention, 
ECG signals are sampled from both a standard 12-lead 
system as shown in FIG. 3 and an orthogonal X, Y, Z lead 
system as shown in FIG. 4. In the Frank lead system de- 
picted in FIG. 4, seven electrodes, designated as H, F, I, E, 
C, A, M are applied to the body at the following locations: 
H, forehead or neck; F, left leg; I, E, C, A, M are located 
along the same transverse level: the fourth intercostal 
space if the patient is supine, or the fifth intercostal space 
if the patient is sitting; I, anterior axillary line; E, center of 
sternum; A, left anterior axillary line; C, at a 45 degree. 



angle between E and A; and M, center of spine. Since one 
can use lead v6 in FIG. 3 for lead A of FIG. 4, and lead v4 
in FIG. 3 for lead C of FIG. 4, there are total 14 electrodes 
needed for the data collection. 

[0034] The acquired ECG data (i.e., high-resolution signals from 
the X, Y Z leads and 12-lead signals) are processed in 
each of a plurality of algorithm modules while data is be- 
ing sampled, as shown in FIG. 5. It is necessary to use 
high-resolution ECG to obtain late potential. Current 
high-resolution ECG analysis already acquires 5-10 min- 
utes of signal, but has not fused the outputs of separate 
algorithm modules as taught herein. 

[0035] | n accordance with one embodiment of the invention, six 
different algorithm modules process respective features In 
the ECG waveform data: late potential, T wave alternans, 
QT dispersion/complexity, QT dynamicity, short-time 
HRV, and intra-QRS notch detection. Each of the six algo- 
rithm modules are available on one and the same plat- 
form. The results of any two or more algorithm modules 
can be fused in combinations that depend on the particu- 
lar clinically relevant group. Diagnosis is based on the re- 
sults of algorithm fusion. 

[0036] The fusion process in accordance with one aspect of the 



invention can be divided into three stages: real-time pro- 
cessing, buffer block processing, and fusion diagnosis, as 
shown in FIG. 6. Each of these stages can, for example, be 
executed by the processor 6 in FIG. 1. The task of the 
real-time processing (step 12 in FIG. 6) is to handle beat- 
by-beat processing needed for all algorithms, which in- 
cludes beat alignment detection for averaging, ectopic 
beat detection for exclusion, on-line signal averaging for 
displaying and on-line parameter calculation for noise 
level monitoring. All of the data, including X, Y, Z Frank 
lead and standard 12-lead ECGs will be stored in a tem- 
porary buffer 14. The data sampling stops in the following 
situations: (a) the noise level of the averaged beat reaches 
a criteria level, e.g., 0.3 uV; (b) a maximum or target 
number of ECG cycles or heart beats to be averaged (set 
by the user) is reached; or (c) the user interrupts the pro- 
cess. The algorithm modules will perform their respective 
processes by using the data in the buffer (step 16 in FIG. 
6). The final diagnosis is based on the fused results from 
all algorithm modules (step 18 in FIG. 6). 
[0037] | n accordance with various aspects of the invention, the 
results of high-resolution ECG waveform data processing 
by two or more algorithm modules can be fused in aid of 



diagnosis. The algorithms can be selected from the group 
of six algorithms described in more detail below. How- 
ever, the broad concept of algorithmic fusion based on 
high-resolution (i.e., signal averaged) ECG data encom- 
passes and anticipates the use of other algorithms in ad- 
dition to those described below. 

[0038] |_ a te PotentialThe QRS late potential module and P wave 
high-resolution module have been implemented in exist- 
ing cardiographs. Substantially the same algorithms, with 
minor modifications, can be a component of the fusion 
diagnostic process disclosed herein. 

[0039] j wave AlternansThe T wave alternans algorithm is based 
on all three leads of the high-resolution ECG data acquisi- 
tion. As shown in FIG. 7, the first step is to divide the 
high-resolution data into even beat and odd beat groups. 
The beats in the two groups are then averaged separately. 
T wave alternans are examined based on three pairs X, Y, 
Z. The algorithm compares intra-variance (i.e., within 
group variance) with extra-variance (i.e., between group 
variance) using a classic statistical analysis procedure 
called ANOVA analysis of variance. In this technique the 
presence of an effect is determined when the variance be- 
tween groups is statistically larger than the variance 



within groups. The intra-variance is the variance of T 
wave morphology of all even beats, or all odd beats for 
the same lead; the extra-variance is the variance of T 
wave morphology between even and odd averaged beats. 
If the extra-variance for two of the three leads are signifi- 
cantly larger than the intra-variance for those leads, T 
wave alternans is positive; otherwise it will be negative. 
[0040] qt DispersionThe known QT dispersion algorithm of Xue 
and Reddy (cited above) is based on 10-sec 12-lead ECGs. 
In accordance with one aspect of the present invention, 
QT dispersion is computed for multiple segments of 
10-sec 12-lead ECGs. The final value of QT dispersion will 
be either average or median values of the QT dispersions 
of all segments. [Since regular QT dispersion is calculated 
from 10-second ECG segments, when 5-10 minutes of 
ECGs are available, we can take either the mean or median 
value of the multiplicity of 10-second-segment QT dis- 
persion values to generate this final value of QT disper- 
sion.] This method provides enhanced reproducibility of 
QT dispersion. 

[0041] qt DynamicityThe QT dynamicity algorithm involves ex- 
amining the beat-by-beat change of the QT interval due 
to changes in heart rate or T morphology. With 200-400 



beats of ECG waveform data from the X, Y, Z leads, the 
processor can compute the short-time QT dynamicity with 
statistically meaningful confidence. 

[0042] short-Time HRVThe heart rate variability algorithm mea- 
sures HRV based on 5 or more minutes of heart beats, and 
is capable of measuring accurately the HF (4-sec rhythm) 
and the LF (10-sec rhythm) components of HRV. An esti- 
mate of the VLF (all rhythms shorter than 5 minutes) com- 
ponent is possible but less accurate. Knowing just the HF 
and LF components is sufficient to give useful information 
about the risk of SCD after myocardial infarction, diabetic 
neuropathy or graft rejection after cardiac transplantation 
among others, as reported by Malik, cited above. 

[0043] intra-QRS AnalysisThe intra-QRS algorithm examines 
small notches inside the QRS complex based on signal- 
averaged high-resolution ECGs. As discussed above, Endt 
has reported that the significant notches found inside the 
QRS complex are correlated with ischemic events before 
and after PTCA. One aspect of the present invention will 
match the notches inside the QRS complex from averaged 
X, Y, Z beats. 

[0044] Fusion DiagnosisThe results of the applying the foregoing 
algorithms to the high-resolution ECG data are then fused 



and used for the purpose of diagnosis. Algorithm fusion 
can be performed in clinically relevant groups. For exam- 
ple, the results from the late potential (either QRS or P 
wave), T wave alternans and QT dynamicity algorithm 
modules can be fused for VT prediction and atrial fibrilla- 
tion. In addition, the results from the QT dispersion and 
intra-QRS algorithm modules can be fused for the pur- 
pose of ischemia detection and prediction. Also the re- 
sults of the Late potential of QRS, HRV, QT dispersion and 
QT dynamicity algorithm modules can be fused for pre- 
dicting sudden cardiac death for post-myocardial infarc- 
tion patients. 

[0045] An algorithm fusion model takes the outputs from differ- 
ent algorithms as its input parameters. This model's out- 
put is the prediction score (ranging from 0 to 100, with 0 
meaning least likely and 100 meaning most likely) for 
particular clinical end points, which can be, for example, 
sudden cardiac death or sustained ventricular tachycardia. 
The predictive model is trained with clinically confirmed 
data for both input and output. The training methods in- 
clude both neural network (nonlinear regression) or linear 
regression. After the model has been trained with known 
outcome output, it can generate the predictive score for 



unknown outcome patients. This fusion model can gener- 
ate better predictive accuracy (with both good sensitivity 
and specificity) than use of a single algorithm output 
without fusion. 

[0046] The advantages of algorithm fusion as disclosed herein 

are manifold. Multiple noninvasive features can be gener- 
ated from the same amount of data collected for high- 
resolution ECGs. The technique provides improved accu- 
racy and efficiency of the noninvasive ECG detection and 
diagnosis processes. Also, algorithm fusion of high- 
resolution ECG data reduces patient/technician study 
time. 

[0047] while the invention has been described with reference to 
preferred embodiments, it will be understood by those 
skilled in the art that various changes may be made and 
equivalents may be substituted for members thereof with- 
out departing from the scope of the invention. In addition, 
many modifications may be made to adapt a particular 
situation to the teachings of the invention without depart- 
ing from the essential scope thereof. Therefore it is in- 
tended that the invention not be limited to the particular 
embodiment disclosed as the best mode contemplated for 
carrying out this invention, but that the invention will in- 



elude all embodiments falling within the scope of the ap- 
pended claims. 



